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Abstract: In case of the discovery of supersymmetry at the LHC, the goal wiU be to 
identify the underlying theory, its fundamental parameters, and the masses of SUSY parti- 
cles. We followed here the approach to reconstruct the decay chains in SUSY events under 
the assumption of common intermediate masses. These masses cannot be extracted from 
each event because of the unmeasured LSP momenta in case of -R-parity conservation. But 
an ensemble of events can be over-constrained, if the decay chains are long enough, such 
that enough mass constraints are available. Here, we present a new method combining a) 
a SUSY mass scan, b) a kinematic fitting based on a genetic algorithm for decay chain 
reconstruction, and c) the usage of angular decay information to suppress the background 
from other SUSY processes. Taking into account the full combinatorial background and 
experimental resolutions in the most difficult case of the fully hadronic decay mode, we 
demonstrate, within one SUSY scenario, that this method can be used to derive a proba- 
bility map of the SUSY parameter space. 
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1. Introduction 



One possible scenario at the Large Hadron Collider (LHC) is the observation of a signifi- 
cant number of events above the Standard Model (SM) background with many jets and/or 
leptons in the final state, as well as large amounts of missing transverse energy ^t- Even if 
the cross sections of these processes are large, as predicted by some models, the identifica- 
tion of the underlying theory is rather challenging. The determination of the masses of the 
involved new particles is not straight forward, if one or more particles escape the detector 
unseen, as it occurs in i?-parity conserving super symmetric (SUSY) extensions of the SM 
Q, where the lightest supersymmetric particle (LSP) is in many scenarios only weakly 
interacting. A further major challenge is the suppression of the combinatorial background, 
which can be of order 10^, if a large number of jets are present. 

There are various approaches to mass determination in SUSY events, which can be 
divided into two major groups: 

(a) Determination of endpoints in invariant mass or other kinematic distributions. These 
observables are typically sensitive to mass differences, and by the combination of a 
few of them, the masses of the involved particles can be accessed |B 



(b) Reconstruction of entire decay chains on event-by-event basis. This is only possible, 
if there are at least as many constraints as unknowns for the problem [E3]-|47]. 



Here, we follow the second approach and, for the first time, take into account the experi- 
mental uncertainty of the measurements for both the event simulation and, by means of a 



kinematic fit, in the event reconstruction |48, 49 1. In order to develop a method which is 
able to deal with high combinatorial background, we demonstrate our method on the fully 
hadronic decay mode of qg production with seven jets in the final state (see Fig. |l|). 

This document is structured as follows: after the introduction, we summarize the idea 
of constrained kinematic fits and present a general technique based on a genetic algo- 
rithm, which simultaneously reconstructs the unmeasured particles and finds the correct 




Figure 1: Feynman diagram of a possible 
fully hadronic SUSY event. 



Figure 2: Feynman diagram of a semilep- 
tonic tt event. 
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jet combination (section ^). The method is tested and vahdated for the reconstruction of 
semileptonic events (section]^. Then, we show the performance of the algorithm for the 
example of one particular SUSY scenario (section The presented technique is finally 
applied to determine the masses (section before we conclude. 

2. Event reconstruction using kinematic fits 

The reconstruction of complete events, including unmeasured particles like neutrinos or 
LSPs, is possible, if there are more constraints than unknowns for the problem. In case 
of cascade- like decays of heavy particles, the invariant masses of the intermediate particles 
and the momentum balance can be used as constraints. As an example, the neutrino 
momentum components (three unknowns) of a semileptonic ti event (see Fig. ^) can be 
reconstructed using the transverse momentum balance (two constraints, T,px = and 
T,py = 0, longitudinal boost unknown), the invariant W masses (two constraints, but 
only the one of the leptonic branch constrains the neutrino), and the invariant top masses 
(two constraints, or only one in case of equal but unknown masses of the top quark). 
In this way it is possible to formulate an over-constrained problem, which can be solved 
by a constrained least square fit, taking into account the experimental uncertainties, and 
therefore not only reconstruct the unmeasured particles, but also improve the resolution 
of the measured ones. 

A common technique to solve such a constrained least square problem is the method 
of Lagrangian multipliers (LM) |^^. Given N measurements yi with the errors fij, the best 
values fi are to be determined, by minimizing the following expression of squared residuals 

s = x' = Y.r-^) (2-1) 

i=l V '^4 / 

subject to M constraints of the form 

Cj(/i, . . . , Jn+p) = for P unmeasured parameters. (2-2) 
The optimal values of fi yield an extremal value of the Lagrange function 

M 

L = S + 2-^Xj-Cj (2.3) 

i=i 

with the auxiliary parameters Xj, called Lagrangian multipliers. The determination of the 
optimal values of / and A are only then analytically possible, if the constraints Cj are 
linear in /j. In particular, this is not the case for invariant mass constraints which are 
quadratic or even more complicated, depending on the representation of the 4-momenta, 
as {px, Py, Pz, E) or {pt, 1], (p, E), where the latter is closely related to the experimental 
measurements and therefore less correlated. Usually, the constraints are linearized, and an 
iterative procedure is applied, to determine the best fitting parameters. Depending on the 
starting values and the number of unknown parameters, it is non-trivial to find a robust 
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algorithm, which improves the solution at each step. Another problem is the definition of 
the convergence criterion, e.g. how well the constraints have to be fulfilled. 

An alternative approach to solve a non-linear constrained least square fit is to interpret 
the constraints cj as additional measurements. The constraints are squared and added to 
the squared residual sum S. If the resulting expression shall be interpreted as an 
function, a proper normalization has to be taken into account. 

In case of an invariant mass constraint of R particles to a mass m with a natural width 
r, the additional term Sm has the form: 

(Ef., Pi)' 



= ^ (2,„r)^ P-"' 

Here, the approximation of normal distributed masses is assumed, to interpret the resulting 
expression as a term. If for different events the masses for intermediate particles at a 
particular position in the chains are not equal, e.g. because of imperfect mass degeneration 
of squarks, the choice of F has to be done carefully (see below). 

In a similar way it is possible to construct a term for the transverse momentum balance 
of the involved particles 



2 

Eall particles i _ n\ 
i=l ^x/y ^ 



St = ^ ^ (2.5) 

The uncertainty of the transverse momentum balance at, assumed to be equal for x and 
y direction, has to be chosen with care. In case of no initial (ISR) and final (FSR) state 
radiation the balance is almost perfect, due to negligible intrinsic transverse momentum 
in the beam. But in realistic scenarios some energy deposition outside the acceptance is 
breaking the balance, and at has to be given a larger value. The choice of F and at has an 
influence on the interpretation of the fit probability. The expression to be minimized is 

S' = S + Sm + St, (2.6) 

which can be done by various algorithms. For reasons which are discussed later, we choose 
an implementation of a genetic algorithm (GA). 

For a given combinatorial hypothesis, the GA searches for the optimal parameter values 
in the following way: All momentum components of all particles (measured and unmea- 
sured) of an event are interpreted as the genome of an individual, such that a complete 
event hypothesis is one individual. The fitness function, which defines if a individual is 
"fitter" than another one, is given by S' in Eq. ( |2.6| ). Starting with a population of indi- 
viduals, which can for example be created by smearing the measurements of the event ^, 
the following algorithm is applied: 

(1) Choose a few fittest individuals, according to their fitness function. 



^The starting values of the unmeasured particles are obtained by assuming the same flight direction as 
a set of measured particles, and choosing the smallest momentum, which fulfills the given mass constraint 
within the experimental resolution. 
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(2) Create a new large generation of individuals by choosing randomly two parents. A 
new individual inherits the momentum components of each final state either from one 
or the other parent. Each momentum component has a small probability to mutate, 
which means it is to be smeared, according to a Gaussian distribution with a given 
width. This width can be given e.g. by the error of the measurement or some other 
carefully chosen value. 

(3) Check the convergence criterion, e.g. no change of the best individual within the last 
few iterations. If convergence is not reached go back to the first step. 

In summary, we traded the problems of the LM technique (linearization and conver- 
gence criterion) for computer resources, which are needed to evaluate the large amount 
of individuals, to find the best fitting values in the high dimensional parameter space. If 
the values of T and at are chosen correctly, the expression S' for the best individual of 
each event will follow a distribution with P — M (number of unmeasured parameter — 
number of constraints) degrees of freedom. 

One of the main challenges of this analysis is the suppression of the combinatorial 
background. The advantage of the genetic algorithm is, that the jet combinatorics can be 
treated simultaneously with the fit of the parameters. To achieve this, the jet combination 
is interpreted as an additional gene. Like the other genes, this one is inherited from either 
one of the two parents. For the jet assignment, the mutation is implemented by (a) simple 
exchange of two jets or (b) exchange of large parts of the two decay branches. 

The algorithm settings, including a time dependency of the mutation, were chosen to 
increase the efficiency to find the correct event interpretation. 



3. Kinematic fit of tt events 

As a proof of principle, we performed the complete event reconstruction of semileptonic ti 
events (Fig. ^). Assuming the top mass is known, this results in six constraints for three 
unknowns, leaving overall three degrees of freedom (see above). The results of the GA are 



compared to a fit with the LM method |51]. The studied tt sample was generated with 



Pythia [£2[. For simplicity, ISR and FSR are switched off in this event generation. The 
true values of the final states were smeared according to typical momentum and angular 
resolutions for jets expected at the LHC experiments: 

- ' > +0.05^ 



pI we 

^2 = f ^] +0.022 



crl={^] +0.02^ (3.1) 



Only events are fitted, in which all four jets as well as the lepton are within the acceptance 
of the detector, assumed here to be pT > 30 GeV and |r/| < 3.0. The event selection 
efficiency for the semileptonic decay mode is 24%. 
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Figure 3: pr resolution of the fitted neu- 
trino in semileptonic ti events. 
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Figure 4: (f> resolution of fitted neutrino in 
semileptonic tt events. 
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In Fig. 1^ and ^ the pT and </> resolutions of the fitted neutrino are shown for the case, 
that the correct assignment of the jets is known. The LM converges for 98% of the events 
for our choice of the convergence criterion {'}2i=i kjl < M ■ 0.1 and AS" < 10^^). The GA 
does not reject any event by construction. 

In Fig. ^ the corresponding fit probabilities are shown for both algorithms. Both 
distributions are almost flat, as expected. Both algorithms show an increase at lowest fit 
probabilities, due to acceptance cuts and the Breit-Wigner distribution of the mass of the 
decaying particle. 

To study the effect of jet combinatorics, 
the assignment of all jets (6 and light quark) 
to the decay chains is left free to the fit. For 
the GA the optimal assignment is obtained 
during the fit process as described above. 
For the LM all combinations are fitted, and 
the one with the smallest probability is 
taken. For both fitting algorithms the se- 
lection of the best combinatorial hypothesis 
will lead to higher fit probabilities, because 
there is a chance for a wrong combination, 
to be better than the right one, whose prob- 
ability is equally distributed. 

In this case the GA finds for 75.6% the 
correct combination as the best hypothesis. 
For the LM fit this efficiency is 72.9%. This 

difference is partially due to the convergence rate of the LM. Furthermore, in case the LM 
selects the correct combination, the GA finds the same hypothesis in 94.6% of the cases. 
Both algorithms lead to the same combinatorial hypothesis in 93.3% of the events. 

For this event topology, for which the number of constraints is considerably larger than 
the number of unknowns, the overall performance of the GA is similar, or slightly better, 
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Figure 5: Fit probabilities for semileptonic ti 
events. Triangles: probability of LM fit for 
the correct jet assignment. Points: same as above, 
but for GA, interpreting S' as x^- 
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compared to the well established LM technique. It is therefore interesting to investigate 
the GA in addition to the LM for the more difficult case of SUSY event reconstruction, 
with more than one unmeasured particle and a factor of ~ 10^ larger jet combinatorics. 



4. SUSY decay chain reconstruction 



For the application of the kinematic fit to supersymmetric cascade decays, we choose on 
purpose the fully hadronic decay mode of a squark-gluino pair (see Fig. ||) as the "signal" 
process, where the combinatorial problem is most severe. The gluino decays into a squark 
and a jet, which will be referred to as "gluino jet". Both squarks then decay further into 
a jet ("squark jet") and either a X2 o^" ^ xt whose masses are in many scenarios almost 
degenerate. The X2 '^^ xf decay into the Xi (LSP) with an additional W, Z, or Higgs 
boson h which also decays hadronically ("boson jets"). In many SUSY scenarios the fully 
hadronic decay modes have the largest branching ratio, and their reconstruction is most 
challenging due to a huge combinatorial problem. Matching the seven jets to the final 
states of the two decay chains, results in 1 260 possible combinations, not counting the 
irrelevant permutations of two jets of the same boson. 

Since two LSPs escape the detector, resulting in six unknowns per event, at least seven 
constraints are needed to be over-constrained, in case of a given SUSY mass hypothesis. 
These are provided by the two components of the transverse momentum balance, the gluino 
mass, the two squark masses, and the two X2/X1 masses. 

For a first reduction of the combinatorial background, boson jet pairs are selected by 
requiring their invariant mass to lie between 50 and 150 GeV, wide enough to account for 
the W , and h masses. The masses of the squarks of the first two generations are assumed 
to be degenerate, as expected in many SUSY scenarios. Only the left-handed squarks are 
considered as signal, since right-handed squarks decay dominantly directly into the LSP. 

The SUSY benchmark point chosen for this analysis has the mSUGRA parameters, 
masses, and properties specified in Tab. [l|. The mass spectrum was calculated with 



SOFTSUSY [53| and the events were generated with Pythia [52]. 



Parameter 


Value 


Particle 


Mass [GeV] 


mo 


230 GeV 




807, 800 




360 GeV 


J/2 
% 


775, 782 


tan (3 


10 


734, 771 


^0 


GeV 


-.1/2 

Qt 


599, 787 


sign^u 


+ 


9 

Al, 2,3,4 
X^2 


851 

144, 271, 475, 490 
273, 487 



Table 1: Parameters and masses of the chosen mSUGRA benchmark point. The total cross section 
is crj^*^ « 1.8 pb for a center of mass energy-^/s = 10 TeV. 



Here, we study both the idealized case without ISR and FSR and the more realistic case 
with ISR and FSR in the event generation (see Tab. ||). Without ISR and FSR, the signal 
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process has exactly seven jets in the final state^. Prom the test sample of 200 000 events, 
corresponding to a integrated luminosity of ~ 50 fb~^, 4.4% remain after the selection cuts 
(no lepton, exactly seven jets with pxj > 30 GeV and \r]j\ < 3.0). The signal to background 
ratio (S/B) of the selected sample is ~ 1/11. In case of ISR and FSR the background is 
much increased, and furthermore only about half of the events have all the jets originating 
from the cascade decays within the selection cuts (see Tab. ^. Note that some background 
events have a rather signal-like decay topology, for example via third generation squarks 
or via heavier neutralinos or charginos. 





without ISR & FSR 
Bg Sig 


with ISR & FSR 
Bg Sig 


selection efficiency 
S/B 
S/B (complete) 


4.2% 29% 
1/10.9 
1/11.3 


6.9% 30% 
1/16.4 
1/33.4 



Table 2: Selection efficiencies and S/B; "complete" means that all jets from the cascade decays 
are within the selection cuts. In not "complete" events at least one jet of the cascade decay is 
outside the acceptance cuts, but instead an ISR, FSR or underlying event jet is selected. 

In a first step, we aim for the complete reconstruction of SUSY events for a fixed set 
of SUSY masses. Therefore, in this section the masses are set to their true values, while in 
section ^ a scan is performed to determine the masses. 

4.1 Kinematic fit 

In the following, we present the fit results including ISR and FSR in the event simulation. 
Due to ISR, the transverse momentum balance of the qg decay products has a width of 
at ~ 40 GeV (see Eq. ^). To compensate for this, the measured jets outside the acceptance 
are added to the momentum balance, improving the above number to at « 15 GeV for the 
fit. For the SUSY mass constraints, the natural width as well as effects from imperfect mass 
degeneration are taken into account (Fg = 10 GeV, Tq = 15 GeV and F-p = 5 GeV). 

In order to compare the performance of the GA and the LM, we study the pT and 
(f) resolution of the fitted neutralinos (see Fig. ^ and ^). The fraction of events for which 
the GA finds the correct jet assignment as the best hypothesis is 42% which has to be 
compared with 1/1260 in the case of random combinatorics. In ~ 10% of the remaining 
events, the best hypothesis is given by the exchange of the decay branches from the two 
squark decays. In this case all mass constraints are also fulfilled, except for the gluino 
mass constraint, which depends on the relativly soft gluino jet^ with its larger relative 
uncertainty (see Eq. |3.lD . In the remaining cases wrong hypotheses are selected due to 
combinations with an accidentally small S' . 

^In our simple detector simulation we neglect the possibility of overlapping jets, which may lead to less 
jets in the final state. The underlying event of multi parton interaction may also lead to additional jets 
within the acceptance, however the number of such events is small (< 4%). 

^Due to the relative small mass gap between the gluino and the squarks, in comparison to the gap 
between squarks and xVxf > the squark jet is on average harder than the gluino jet. 
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Figure 6: pr resolution of fitted neutrali- 
nos. Only events are shown where the best 
hypothesis has the correct jet assignment. 
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Figure 7: resolution of fitted neutralinos. 
Only events are shown where the best hy- 
pothesis has the correct jet assignment. 



The definition of a convergence criterion for the LM is difficult. Choosing a setting 
which yields a small bias and resolutions in pT and similar to the GA, we obtain the 
same "purity" (41%) of correct jet assignments, however at the cost of a much reduced 
convergence rate of ~ 50%. Due to its higher rate, the GA is the algorithm of choice for 
the further studies. 

In Fig. ^ the fit probability p for the 
GA is shown. Above p > 0.1 the distribu- 
tion looks reasonably flat for the signal, as 
expected. Possible differences to a flat dis- 
tribution are the jet combinatorics, which 
leads to a shift towards higher probabili- 
ties, and under- or overestimation of un- 
certainties. In addition, the distributions 
are shown for the total background, as well 
as for "signal-like" background. The latter 
consists of either signal events where one 
jet of the cascades lies outside the accep- 
tance cuts and instead an ISR or FSR jet 
is found, or decays with the correct decay 
topology but involving squarks of the third generation. The number of such events is of 
the same order as the number of signal events. They also show a rather flat probability 
distribution for values beyond 0.1. In contrast, the background events show a clear de- 
viation from the uniform distribution with a strong increase towards lower values. As a 
consequence, a cut on p can be used to improve the signal to background ratio. In this 
example, the total S/B of 1/33 is improved to 1/11 (1/8) for a cut at 0.1 (0.3) at the cost 
of reducing the signal efficiency by 79% (68%). It is noticed, that beside the rather simple 
selection cuts no further reduction of the SUSY background has been attempted here. 
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Figure 8: Fit probability for best event hypothe- 
ses using the GA. 
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Figure 9: Angular distribution of the bo- 
son with respect to the flight direction of the 
/ X2 ill the rest frame of the decay, for the 
GA fit. 
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Figure 10: Angular distribution of the 
squark jet with respect to the flight direc- 
tion of the squark in the rest frame of the 
decay, for the GA fit. 



4.2 Angular distributions 

So far, we have only utilized invariant mass and transverse momentum balance constraints. 
Further information of the event can be used, like for example decay angles. In Fig. |9| and 
10 the cos Q* distribution of the decay products in the rest frame of the decaying ^ 



and squark are shown, where Q* is the fitted angle with respect to the flight direction of 
the decaying particle. In case of a spin-0 or unpolarized particle, the expected distribution 
is uniform, as can be seen for the correct jet combinations. The acceptance cuts on the jets 
influence mostly the distribution for the xf /X2 decay: The lower cut prefers bosons 
decaying into the flight direction of the SUSY particle, since then the two jets have on 
average higher p-p. This leads to a slight increase in the cos0* distribution towards -|-1. 
For the background consisting of other decay modes the distribution is very different. 

The normalized distribution of the background is parametrized by a polynomial and 
then used as a probability density / for a likelihood ratio. The signal distribution is 
assumed to be flat. The probability density to be signal-like of an event hypothesis with a 
set of N decay angles cos Q* can then be expressed by: 

ro_ 1 TT /sig(cosgr) . 

Cnorm fj /sig(cOS 0*) + /bg(cOS 0*) ^ ' ^ 

The product L 

L = p-LR (4.2) 

is used to define the best hypothesis, including angular information of the event. In this 
case L helps for a further suppression of the SUSY and SM background. In our example 
a cut on L with a signal efficiency of 77% leads to S/B of 1/9, while a cut on p at 
0.1 with a similar efficiency of 78% yields 1/11. For the reduction of the combinatorial 
background, the performance of L and p are similar. It should be pointed out, that the 
probability distribution of the background depends on the model, which was used for this 
study. Further studies have to be made to reduce this model dependence. 
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5. Constraining SUSY masses 



In the following the mass reconstruction is studied using the GA including full combi- 
natorics but ignoring other SUSY processes. For this a scan of m{q) and m{xf/x2) is 
performed within the allowed kinematic range, fixing 171(0) and 'm{xi) to their true values. 
Although two masses cannot be constrained by a single event, a set of events leads to an 
over constrained system. The true mass hypothesis is expected to result on average in a 
better fit probability. There are various measures to quantify this behavior: 

• Number of converging events iVconverged- In the case of the GA "converging" means 
that the fit probability is above a given threshold Pcut- In the following we choose 
Pcut = 0.01. 



Average fit probability of converging events: This measure has the disadvantage that 
single events with a good fit probability for a wrong mass hypothesis can suggest 
"good" regions in the mass space, although only a few events are converging at all. 

Fractional event count ATf^g^: This is defined as the sum of the fit probabilities of 
all events. This quantity has only a small dependence on the number of converged 
events, since not converged events do not contribute due to their small probability. 
However, a statistical interpretation is difficult. 

Normalized likelihood of all events: The product V of the fit probabilities pi for all 
events is a natural definition of the overall likelihood of a given mass hypothesis. Un- 
fortunately, the number of converging events is not constant for all mass hypotheses. 
The likelihood can be normalized to the number of events as follows: 



^rr,i-wi-'rcrnr\ \ ^Converged / ^ -^converged 



n = exp — logpi (5.1) 

i=l / V^Vconverged -^^ 



In regions where the number of converging events is small, T^norm 

suffers from the 

same fluctuations as described for the average fit probability. 

• Regularized likelihood: To avoid the problems of normalized quantities, each not 
converged event can be assigned a small probability Pcut- In this case, the same 
number of events are used for each mass hypothesis, and V does not have to be 
normalized. The logarithm of V has numerical advantages and can be interpreted 
statistically. 

log = XI ^^iPi,Pcut) (5.2) 
i=l 

Note, that the value of Pcut is not critical since the statistical interpretation depends 
on A log V relative to the maximum, and close to the maximum variations are negli- 
gible. 
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Figure 11: A'fi.ac of events for which the GA 
finds the correct jet assignment. fSR and 
FSR are included. 
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Figure 12: A^frac of events for which the GA 
finds the wrong jet assignment. ISR and FSR 
are included. 



In Fig. (12) the fractional event count A'^f^ac is shown in the xf /X2 ~ Q mass plane 
for the events where the fit finds the correct (wrong) jet assignment. For the events where 
the correct hypothesis was found, the true masses lie within 10 GeV around the maximum, 
in contrast to the wrong combinations which prefer a larger xf/x2 ™ciss. We attribute 
this to the larger combinatorics for the xf /X2 heavier mass hypotheses, as well as to jet 
threshold effects at low mass hypotheses. In principle, if all the uncertainties are known, 
this bias can be extracted from Monte-Carlo simulations. 

Fig. |l^ shows A log V without combi- 
natorial background. As for the case of 
Afrac the true masses are found with a pre- 
cision similar to the width F of the mass 
constraints. If the combinatorial background 
is included this picture changes as shown 
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in Fig. |lj. A clear correspondence between 
the maximum of A log V and the true squark 
mass of 800 GeV is visible. However, there 
is a systematic shift towards larger xf/X2 
masses, with respect to the true value of 
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Figure 13: A log 7^ of signal events, including 
ISR and FSR, and without combinatorial back- 
ground. 



270 GeV. In Fig. 15 AlogP is shown for a 
sample without ISR and FSR. In this case 
the momentum balance is a much stronger 
constraint, since additional jets outside the 
acceptance cuts are not present. As a consequence the GA finds in 48% of the events the 
correct jet assignment, compared to 42% with ISR and FSR included, and the systematic 
shift towards larger Xi /X2 masses is significantly reduced. 

If the SUSY background is included, the interpretation will be more difficult, depending 
on the signal to background ratio. A combination with other observables, like event rates 
or mass edges from other final states, is statistically possible. 
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Figure 14: AlogT' of signal events. The 
true masses are m{x^/x2) ~ 270 GeV and 
m{q) « 800 GeV. ISR and FSR are in- 
cluded. 
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Figure 15: Same quantity as on the left 
hand side, but no ISR and FSR present. 



6. Conclusion 

We have presented a novel approach to perform constrained kinematic fits based on a 
genetic algorithm, and applied it to the complete reconstruction of supersymmetric decay 
chains of qg production with two unmeasured particles. The fully hadronic decay mode 
with seven jets in the final state was chosen, since it has in many SUSY scenarios a large 
branching fraction, and the large jet multiplicity leads to a huge combinatorics, which is a 
viable test for the developed fitting technique. In comparison with a conventional kinematic 
fitting based on the method of Lagrangian multipliers, we observe the genetic algorithm to 
converge for a factor of two more events and to find in ~ 45% of the events the correct jet 
assignment. In addition, decay angles were used to improve the rejection of other SUSY 
processes. In this way more information than only invariant mass constraints were used. 

A SUSY mass scan was performed, and a regularized likelihood was used to derive a 
probability map of the SUSY mass space. We note that squark masses are reconstructed 
quite well, if samples with only signal events are used. The xf /X2 masses are reconstructed 
with a systematic shift towards larger values. This effect is mainly due to events where the 
fit finds a wrong jet assignment, but can be corrected using Monte-Carlo simulations. 

For samples including the full SUSY background further selection cuts to improve the 
signal to background ratio have to be studied. The application of this method to leptonic 
decay modes will be of particular interest, due to reduced combinatorics, less SUSY and 
SM backgrounds, as well as more accurate final state measurements. 
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